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ABSTRACT 

We present Herschel-SPIRE observations at 250-500 yum of the giant elliptical galaxy M86 and examine the distribution of the resolved cold dust 
emission and its relation with other galactic tracers. The SPIRE images reveal three dust components: emission from the central region; a dust 
lane extending north-south; and a bright emission feature 10 kpc to the south-east. We estimate that ~ 10^ M© of dust is spatially coincident with 
atomic and ionized hydrogen, originating from stripped material from the nearby spiral NGC 4438 due to recent tidal interactions with M86. The 
gas-to-dust ratio of the cold gas component ranges from ~ 20 - 80. We discuss the different heating mechanisms for the dust features. 
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1. Introduction 

Studies of cold dust in elliptical galaxies has been limited to 
date by the lack of high-resolution, long wavelength spectral 
coverage. In particular, the origin of far-infrared (FIR) emission 
in these systems is still a controversial issue, with evidence of 
dusty disks favouring a stellar origin (e.g. Knapp et al. 1989) 
and other systems originating from mergers with dust-rich galax- 
ies (Leeuw et al. 2008). The unprecedented resolution and sen- 
sitivity of the recently launched Herschel Space Observatory 
(Pilbratt et al. 2010) allows us to address long-standing issues 
such as the origin and quantity of dust in ellipticals. 

One of the most well-known IR-bright ellipticals is the giant 
Virgo Cluster member, M86, at a distance of 17 Mpc (Mei et al. 
2007). Two dust features were detected with IRAS, one coinci- 
dent with the galaxy and another to the north-west, thought to 
be coincident with an X-ray plume of gas (Knapp et al. 1989; 
White et al. 1991) and originally attributed to dust stripped from 
M86 due to its motion through the cluster. Higher resolution 
data from ISO revealed two dust peaks within M86 suggest- 
ing a massive dust complex (Stickel et al. 2003, hereafter SOS). 
They proposed a tidal origin, also supported by absorption fea- 
tures attributed to dust stripped from the nearby dwarf galaxy, 
VCC 882 (Elmegreen et al. 2000). The discovery of atomic gas 
offset from the centre of M86 and decoupled from its stellar disk 
supports the tidal scenario (Li & van Gorkom 2001). More re- 
cently, Kenney et al. (2008, hereafter K08) detected strong Ha 
features extending from the nearby spiral NGC 4438 at 23' away. 



* Herschel is an ESA space observatory with science instruments 
provided by European-led Principal Investigator consortia and with im- 
portant participation from NASA. 



to within 1' of the nucleus of M86. The distribution and veloc- 
ity of the ionized gas provides clear evidence for tidal interaction 
between these two giants (K08, Fig. 1 Cortese et al. 2010). In this 
scenario, we are observing debris from the spiral left in the wake 
of the collision with M86 with ~ lO^M© of cold gas removed 
from NGC 4438 (K08). The stripped material is then heated by 
the hot interstellar medium (ISM) of M86 or shock fronts from 
the interaction. Although tidal stripping from NGC 4438 is sup- 
ported by the atomic and ionized gas distribution, the origin of 
the dust responsible for the FIR emission and its heating mech- 
anism is unclear. Here we present submm observations of M86 
with Herschel-SPTRE (Griffin et al. 2010). Companion observa- 
tions of NGC 4438, are presented in Cortese et al. (2010). 



2. Observations and data reduction 

M86 was observed with SPIRE at 250, 350 and 500 jum during 
Herschel 's Science Demonstration Phase as part of the Herschel 
Reference Survey (Boselli et al. 2010a). Eight pairs of cross- 
linked observations were taken in scan-map mode with scan- 
ning rate 30"/sec. The data were processed following the de- 
tailed description given in Pohlen et al. (2010) and Bendo et 
al. (2010a). Calibration methods and accuracies are outlined in 
Swinyard et al. (2010). The measured 1 cr noise level is 5, 6 and 
7 mJy beam"' at 250, 350 and 500 yum with beam size 18, 25 and 
37"; the noise is dominated by background source confusion. 



3. Results 

The three-colour SPIRE image is shown in Fig.[T]with dust fea- 
tures labeled (following the terminology of S03) along with an 
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Fig. 1. Le/f.- 7? -band image of M86 region with Ha (red) and X-ray (green) contours. Right: Three-colour image at 250, 350 and 
500 ;um, smoothed to 500 yum beam (as indicated by dashed circle in upper right of the figure). The ellipse indicates the optical halo 
ofM86(8.9'x5.6'). 



optical image of the same region with X-ray and Ha contours. 
In the SPIRE image, there are a number of unresolved sources, 
and at least five extended features. Careful comparison of large 
scale FUV and IRAS maps of this region show that the cir- 
rus emission is extremely low ruling out a Galactic origin. In 
Fig. 12] we focus on the central 13' x 8' region and compare with 
structures seen at other wavelengths. Towards the south, M86- 
S is coincident with a number of clustered 24/im sources. In 
the north-west, the bright source M86-NW (originally associated 
with M86's X-ray plume) has no optical or UV counterpart and 
peaks at A < lOOjum. The bright feature M86-F1R5 is coinci- 
dent with an optical galaxy (VPC 463) and, like M86-F1R4, has 
a similar flux ratio to M86-NW. Extending north from M86, a 
faint (2-3cr) filament-like submm structure appears to be coinci- 
dent with blueshifted Ha emission (Fig.|2li attributed to ionized 
debris from the incoming trajectory of the NGC 4438 collision 
(Trinchieri & di Serego Alighieri 1991; Finoguenov et al. 2004; 
K08). However, this submm filament is also coincident with a 
number of distant 24 /im sources. We therefore suggest that this 
feature along with M86-S, M86-NW, M86-F1R5 and M86-FIR4 
originate from background sources unrelated to M86. 

The submm morphology within the inner 2' of M86 is sur- 
prisingly complex and differs from the smooth distributions seen 
in the 24 /im, UV and optical images. We see four distinct fea- 
tures in the SPIRE maps labeled as M86-N, M86-off^, M86-SE 
and 'lane' (Fig. [1). M86-N is a bright unresolved source which 
is included in the photometry due to its coincidence with a Ha 
knot in this regioifl The feature M86-SE (7.7 x 2.3 kpc) is the 
brightest feature within the optical halo and is ofi'set from the 
centre by 1.9' (lOkpc). To the immediate north of this feature 
extends the bright iane' structure. These features are not seen in 
the Spitzer 24 or 70jum data (although the latter exhibits severe 
striping effects) but are detected with ISO at 90, 135 and 180/im 
and Spitzer 160//m (Fig.|2]). Although there are four 24fim point 



' We note that M86-N has similar features to other point sources in 
the crowded field and could be a background object. It contributes 15% 
of the flux measured in the aperture for M86-off. 



sources coincident with the M86-SE structure, these are likely 
to be background sources since corresponding features are not 
seen in the same position in Ha nor are they coincident with 
the peaks in the submm emission. These four features (Fig. |2]) 
are spatially coincident with peaks in the redshifted ionized gas 
(M86-SE: -120km s"', M86-ofr: -200- 350 km s"') extending 
from NGC 4438 (K08) and with atomic Hi (Li & van Gorkom 
2001). The submm and Hi peak in M86-SE is coincident with 
an Ha 'hole', but this is likely due to the lower resolution of 
the former. M86-SE is also contained within an X-ray bound- 
ary tracing the X-ray halo (Finoguenov et al. 2004; Randall et 
al. 2008). The strong spatial correlation between the cool atomic 
gas, cold dust, ionised gas and hot X-ray boundary shown in 
Fig.|2]seems to suggest that the dust originates from gas stripped 
from NGC 4438 and is immersed in the X-Ray halo of M86. 

At first glance, it appears that the 'lane' feature and M86-SE 
are coincident with dusty filaments seen in absorption (A and B 
in Fig. 121 Elmegreen et al. 2000). Careful comparison with the 
SPIRE images show that their B feature is 1 .3' south of M86-SE, 
and filament A is offset by 0.4'. Faint submm emission is seen at 
the southern tip of A and the atomic gas encompasses feature B. 
We are not able to resolve A (at only 6" across in the optical), 
but surprisingly we do not detect significant emission associated 
with B. This could be a result of the absorption features arising 
from foreground dust lanes with low column density, or the dust 
is simply too cold to be detected with Herschel. 

3.1. Dust mass and heating witiiin the optical halo 

We performed aperture photometry on the IR and submm im- 
ages from ISO, Spitzer and SPIRE. The datasets were wcs- 
aligned, and smoothed to the 160yum beam using the appropri- 
ate kernels (Bendo et al. 2010b) and assuming gaussian beam 
profiles for SPIRE. Fluxes were measured using two elliptical 
apertures (Fig. ^ with semi-major, semi-minor axes and posi- 
tion angles of 1.4' x 1.0', 127.5° (encompassing M86-SE) and 
1.0' X 1.6', 31.5° (encompassing M86-off and M86-N) respec- 
tively. We used an empirical 160yum PSF to determine the aper- 
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Fig. 2. Multiwavelength comparison of M86. Top: UV (GALEX) and optical (SDSS) with absorption features A & B. Middle: 
Spitzer 24-160//m and SPIRE 250-500 //m. White ellipses indicate photometry apertures. Bottom: Gas emission shown in negative 
greyscale. Hi with white contours 10 + 25mJybeam ' (courtesy J. van Gorkom). Ho- (courtesy J. Kenney) and X-ray (ROSAT). 
250 /im shown by red contours 9-1-6 mJy beam '. 



ture correction (Young, Bendo & Lucero 2009). The spectral en- 
ergy distributions (SEDs) of M86-ofr and M86-SE are shown 
in Fig. [3] (see also Boselli et al. 2010b). Blackbody functions 
modified with a A^^ emissivity law are plotted for comparison 
with temperatures 15 K and 44 K. Assuming a mass-absorption 
coefficient, /C500 = 0.1 kg"' m^ (Draine 2003), we can estimate 
a rough dust mass from the total 500 jum flux, 5 500 assuming 
Md = S yD^ I KyB(v,T) where D is the distance. The total dust 
mass within the inner 2' of M86 ranges from ~ 2-5x10^ M© for 
temperatures, T^ - \5 - 20 K. Comparing the dust masses with 
the atomic gas mass, we estimate a range of gas-to-dust ratios of 
gid ~ 20 - 80 depending on the temperature. Non-detection of 



CO gas within M86 rules out a significant molecular component 
with < 10''** Mo (Braine, Henkel & Wiklind 1997). The ionized 
gas could contribute a further 10^ M©, but this is difficult to esti- 
mate without knowing the gas densities (K08). These ratios are 
higher than expected in elliptical galaxies (e.g. Fich & Hodge 
1993) and similar to the lower range of values estimated for the 
dust feature extending out from NGC 4438 (Cortese et al. 2010). 

It is interesting to ask what is responsible for heating the 
dust and creating the submm emission seen here. Possible mech- 
anisms include AGN, the interstellar radiation field (ISRF), tidal 
heating and/or the hot X-ray halo. Following Thomas et al. 
(2002), the luminosity L/, required to heat a cloud of dust with 
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grain size a, total mass Ma, temperature T^ and Planck-averaged 
absorption coefficient < Q(a, 7^) >, is given by Eq.[T] 



Lh — 



Ana-' 



< Q(a, Td) > T^crMd 



(1) 



For a < 0.1 /^m at 20 K, we require L/, < 10^ L©. M86 is ra- 
dio quiet and we can rule out significant heating from the AGN. 
Although the dust is not spatially coincident with the optical/UV, 
the ISRF could still be responsible for heating the dust with L/, 
comparable to the B-Band luminosity (Mei et al. 2007). To in- 
vestigate this, we used a radiative transfer implementation of 
the Monte Carlo code skirt (Baes & Dejonghe 2002), which 
models the absorption, scattering and thermal emission of cir- 
cumstellar discs and dusty galaxies. M86 was represented by a 
flattened Sersic model based on parameters from the literature 
(Caon, CapaccioH & D'Onofrio 1993; Graham & Colless 1997; 
Gavazzi et al. 2005) and we used the global SED from Boselli 
et al. (2010b) for the intrinsic model. During the Monte Carlo 
simulation, the mean intensity of the radiation field is calculated 
from UV-MIR wavelengths at every position in the galaxy. From 
this mean intensity, the equilibrium dust temperature, Teq of each 
species of dust grains is calculated using energy balance. We 
predict that Tgq for silicate (graphite) grains at a distance of 1 .9' 
from the nucleus is 13 (19)K; even though the submm peaks 
do not correlate with optical/UV emission, the ISRF of M86 is 
sufficient to produce the submm emission detected here. 

The thermal energy provided by M86's hot. X-ray emitting 
ISM could also provide a significant heating contribution by col- 
lisionally heating the dust, with grains of a < 0.1 fim reaching 
Teq < 18 K (e.g. Dwek 1987). The tidally-heated component 
as traced by Ha also cannot be ruled out as it provides compa- 
rable energy to M86's thermal reservoir over the last 100 Myr 
(~ 10'^ Lo, K08). Thus, the ISRF, X-ray halo and tidal heat- 
ing could all be contributing to heating the dust and it is diffi- 
cult to separate these processes using the UV-submm SED only. 
However, the high resolution images presented here suggest that 
submm emission only originates from regions where we detect 
atomic and ionised gas. Submm emission is not detected from 
regions with just atomic Hi and X-ray emission, suggesting that 
the X-ray halo alone is not responsible for heating the dust. The 
spatial distribution of these tracers therefore favours a scenario 
in which the submm emission originates from dust mixed with 
stripped atomic material and is heated by the tidal interaction. 

In summary, we present submm images of M86 which reveal 
a complex dust morphology with emission detected lOkpc from 
the centre, towards the southeast. The unprecedented resolution 
of SPIRE has revealed, for the first time, a strong spatial corre- 
lation between the cold dust (~ lO*" Mq, g/d of 20 - 80) and the 
warm ionized gas in M86. This result strongly favours a scenario 
whereby the dust in M86 originates from material stripped from 
the nearby spiral NGC 4438. We investigate the different heat- 
ing mechanisms responsible for the dust emission detected by 
Herschel. Intriguingly, although we cannot rule out the stellar ra- 
diation field of M86, the strong correlation between submm and 
Ha emission suggests the cold dust is heated by the same mech- 
anisms responsible for ionizing the gas stripped from NGC4438. 
If so, tidal heating is likely to be responsible for the dust emis- 
sion. Further modeling is required to provide a definite answer 
on the origin of the submm features revealed by SPIRE in M86. 
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Fig. 3. The SEDs for top: M86-SE and bottom: M86-off. Solid 
lines are modified blackbody functions with cold temperature 
15 K and 44 K. Spitzer fluxes are blue triangles, ISO (stars) and 
SPIRE (circles). IRAS data (green triangles) are measured in a 
single aperture of radius 136". 
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